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We study the consistency of two Higgs doublet models in light of the new bosonic particle discovery 
at the LHC. We work within a general setup that we call the 2HDM-X, in which the quarks couple 
to both scalar doublets with aligned couplings such that flavor-changing neutral currents are absent 
at tree level. The 2HDM-X encompasses the traditional Type 1, Type II, lepton specific, and flipped 
models, but also provides for more general possibilities. The best fit to the current data with 
a suppressed rr signal and a bb signal of Standard Model strength is given by the 2HDM-X with 
specific parameter choices; however, very good fits are also obtained within the lepton-specific model 
and a democratic model, the 2HDM-D, in which both the up- type and down-type quarks couple to 
each doublet with equal strengths. The approach provides a general framework in which to interpret 
future LHC Higgs data within extensions of the Standard Model with two Higgs doublets. 

PACS numbers: 12.60.Fr, 14.80.Ec 



Introduction. The recent discovery of a Higgs-likc par- 
ticle with a mass of order 125 GeV at the Large Hadron 
Collider (LHC) may be the most significant crowning 
achievement in particle physics in several decades. If 
this particle is indeed a Higgs boson, it is of paramount 
importance in confirming spontaneous electroweak sym- 
metry breaking as the origin of mass. The properties of 
the Higgs boson at colliders are also sensitive to physics 
beyond the Standard Model (SM). Knowing or exclud- 
ing additional new particles that affect the Higgs proper- 
ties can increase our understanding of the nature of elec- 
troweak symmetry breaking. One of the best-motivated 
extensions of the SM is to add an additional Higgs dou- 
blet, resulting in two Higgs doublet models (2HDM's) 
(see for a recent review). Specific realizations of 
2HDM's have already been studied in light of the Higgs- 
like particle discovery [H-Q. 

There are several well-known 2HDM's that correspond 
to different ways of coupling the Higgs doublets to the 
SM fermions. The models can be classified into two 
general categories depending on whether they result in 
Higgs-mediated flavor-changing neutral currents (FCNC) 
at tree level. Within the first class of models, tree-level 
FCNC are forbidden by global symmetries such as Peccei- 
Quinn symmetry Q or discrete symmetries. This class of 
models include the Type I (2HDM-I), Type II (2HDM- 
II), lepton-specific (2HDM-L) and flipped 2HDM, which 
can be distinguished based on Yukawa coupling measure- 
ments 0. The second class of models, which are com- 
monly known as Type HI models, include Yukawa cou- 
plings of the up- type and down-type quarks to both Higgs 
doublets. This class of models generically have tree-level 
FCNC and are highly constrained by flavor observables 
which are consistent with the SM. 

Motivated by the experimental hints of deviations of 
the Higgs boson properties from the SM predictions 
and the agreement of the flavor observables with the 



SM. we explore a simple 2HDM model framework that 
captures the dominant effects of Type III models in 
modifying the Higgs properties and at the same time has 
no tree-level FCNC. In this setup, which we denote as 
the 2HDM-X, the quarks couple to both Higgs doublets 
with a flavor structure that is governed by the same 
Yukawa matrix This coupling choice guarantees 

the absence of tree-level FCNC due to fine-tuning, and 
thus a priori FCNC may reappear at higher-loop level 
with appreciable strength. However, we choose this 
general framework because our principal guidance is 
the experimental Higgs data, which points to specific 
features that are not yet understood. Before proceeding, 
we also comment that the label 2HDM-X we use here 
should not be confused with previous literature which 
also use this title in other contexts, such as [§]. 

The 2HDM-X. In this setup, there are two complex 
Higgs doublets $d = ($°, $^) with hypercharge Y = —1 
and = {^t,^u) "^ith Y = +1. After electroweak 
symmetry breaking and assuming that the tree-level 
Higgs potential conserves CP, we have 

*d = [ {vd + + 1 /V2, $-] , 

= [$+, (iv + </);+*C)/V2] , (1) 

in which and are the usual electroweak vacuum 
expectation values (VEV's), with = w„ + v'^ — 
(246 GeV)2. The ratio of the two VEV's is defined as 
usual to be tan /3 = Vu/vd- In the scalar spectrum, there 
are two neutral scalars h and H, one pseudo-scalar A, 
and one charged scalar H^. Here we mainly consider 
the phenomenology of the lightest CP-even neutral scalar 
h, but we will also include the effects of the charged 
scalar on the h properties as well as the experimental 
constraints on H, A and H^. The lightest CP-even 
scalar h is given as usual by ft, = sin a + (f>^cosa 
with the mixing angle —tt/2 < a < tt/2. In the de- 



2 



coupling limit with ruh <C ~ mu ~ fTiH±, one has 
cos {(] - a) = 0{v'^/m\) (T^]. The tree- level couplings of 
h to the W and Z are 



= 9v m-v sin (/? - a) 



(2) 



in which gv = 2my/u for V — W or Z . The Yukawa 
couplings of the Higgs doublets to the quarks and the 
charged leptons are given in the 2HDM-X without loss 
of generality as follows: 

- C = Vu UR (cos 7„ + sin 7„ Ql 
+ yddB. (cos7d *d + sin7d *„) Ql 
+ yteR^dLL + h.c. , (3) 

in which yu,d,t are 3x3 Yukawa matrices and ^u,d = 
itT2^.u d' ^'^'^ have neglected the neutrino interactions 
with the Higgs for simplicity. In the above, we have used 
the freedom to redefine the two linear combinations of 
$u and [m to eliminate the coupling of the leptons 
to $11. Hence, in this basis there are two mixing angles 
^u.d for the up-type and down-type quark couplings. 

The 2HDM-X parameter space reduces in specific lim- 
its to several familiar 2HDM's, as shown in Table HI It is 
clear that the 2HDM-X will allow us to cover a broader 
range in the Higgs sector than any of these familiar limits. 
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TABLE I: The values of mixing angles 7u,d for different types 
of 2HDM's. "L" means the lepton-specific 2HDM (2HDM-L). 
The "flipped" 2HDM (2HDM-F) is similar to the 2HDM-II 
except that the charged leptons couple to the same Higgs as 
the up-type quarks. The "democratic" (2HDM-D) refers to 
the case in which both doublets couple equally to up and 
down quarks. The last column contains the best fit values of 
the 2HDM-X in light of the current LHC data. 

The lightest CP-even neutral scalar couplings to 
fcrmions, in units of the SM coupling, take the form 
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in which we have shown only the dominant tree-level cou- 
plings to the third-generation fcrmions. For the charged 
Higgs i?^, the Yukawa couplings are given by 



9H-T+V 



V2 

V 

V2 



9H-tb = — [mt cot {13 + ^u)PR + mbia.n {13 -jd)PL] 

V 



m-r tan/3 Pi\ 



(6) 



Finally, the scalar potential is key for determining the 
IiH^H" coupling that arises in the /i — >■ 77 decay. 
The scalar potential parameters can be traded for the 
masses and angles Mh, Mh, M^, M[j± and a, such that 
the charged Higgs coupling can be written as 



9hH+H- 



= 2 



cos (a + 13) Ml - 
sin 2/3 V 



sin(/3 — a) 



Ml 



(7) 



where = 2m 12 / sin 2 f3 is the common mass scale of 
the heavy Higgs sector, and mfj is the mass-term of the 
*I'i*I'2 + h.c. term of the Lagrangian. We note here that 
our derivation of this coupling has the opposite sign as 
compared to that of [ld| . 

From Eq. ^ and Eq. (O, we see that the SM Higgs 
boson couplings can be recovered by choosing /3 — a+TT/2 
and 7u = 7rf = 0. Gcnerically, the couplings to weak 
gauge bosons are smaller in the 2HDM-X than the SM 
values. However, the couplings to fcrmions can be cither 
larger or smaller than the corresponding SM couplings. 
Also, unlike the discrete 2HDM's, the couplings to the 
t-quark, the 6-quark and the r-lepton are independent of 
each other. For example, by choosing a = and non- 
zero values of 7„ and 7^, one can even have a vanishing 
h ^ TT partial decay width. As we will see later, it is 
this freedom that provides a better fit for the 2HDM-X 
than other 2HDM's. 

In order to compare with the experimental measure- 
ments, it is necessary to determine the branching frac- 
tions of various Higgs decay channels in this setup. First, 
we determine the total width of h, which is given by 
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cos2 (/3 - 7d) 

where for illustrative purposes we have assumed that the 
total width is always dominated by the bb and W~^W~ 
channels and we have neglected subdominant contribu- 
tions from TT and gg. In our numerical calculations, 
we will include all the subdominant channels in the to- 
tal width calculation. Given the parameter space free- 
dom, the total width can be much larger or smaller 
than the total width of the SM Higgs boson. The ra- 
tio of branching ratios for several important channels. 
Pi = BF,(2HDMx)/BFi(SM), then take the form 
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(11) 



47 sin (/3 + 7„) 

we show here the result in the limit to;, 
2to4, 2Mh/, but in our numerical calculations we use the 
full one-loop formula including the contribution. 
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FIG. 1: Ranges of a and tan /3 consistent with the experimen- 
tal data. A random selection of points in our Bayesian fit are 
presented as red dots. The point with the best overall fit is 
marked by the blue x. Data consisting of decays to 77, VV , 
bb and rr constrain the parameter space to regions of black, 
blue, red, and purple colors, respectively. 



The experimental data are typically presented as 
ratios of the production cross section times branching 
fraction divided by the corresponding SM value for 
the production type "i" and the decay channel ".7"). For 
the Higgs resonance and using the Breit-Wigner formula, 
the can be related to the pi by = Pj Pi for 

the exclusive channel. For the inclusive channel, we sum 
over the gg, VV and Vh production modes. 



Fit to the Higgs data. Wc perform a Bayesian fit to 
the available LHC and Tevatron data folloviring the same 
method as [13) ■ We require all scalar potential parame- 
ters to be perturbative with A? < Att. From the present 
experimental measurements 13-1^ (see also 0)i 8- sim- 
ple combination of inclusive and selected exclusive chan- 
nels provides the following distilled measurements: 
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(12) 

in which pp, gg, VV and Vh indicate inclusive, gluon fu- 
sion, vector boson fusion and associated production of the 
Higgs boson, respectively. Furthermore, we constrain the 
H state with the excluded signal rate over a mass range 
up to 500 GeV 

It is well known that a light charged Higgs boson can 
increase the ft, — 77 rate. Therefore, wc include the 
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FIG. 2: Posterior probability of the heavy Higgs masses. To 
achieve the enhanced /i — > 77 rate, the charged Higgs is typ- 
ically light. Consistency with electroweak precision observ- 
ables forces the CP-odd to be of similar mass. 



values of the precision electroweak observables S and T 
and the flavor changing decays B^^Xs + ^ + Xat next 
to leading order in our fit 0. We require S = 0.04±0.09, 
T = 0.07±0.08 with an 88% positive correlatio n 1 171 and 
BF(Bo Xs+i + X)^ (3.55 ± 0.26) x 10"^ 

In addition to the enhancement in the 77 channel, a 
curious feature of the measurements is an apparent sup- 
pression in the rr channel. We will sec that this can be 
accommodated in the traditional lepton-specific model as 
well as the democratic and general models. 

In Fig. [U we show the bands in a and tan/3 near 
the best fit point associated with the collider measure- 
ments in Eq. ((T2|) . The solid lines indicate the central 
value while the shaded region indicates the Icr band of 
the measurements combined over all production channels 
[ii{VV) = /i(77) = 1-9 



1.3tJ;^ and 



p{tt) = 0.0tQ;o]- The Bayesian fit is represented by the 
red points. Therefore, the measured VV , 77 and rr rates 
provide a good fit for large tan (i and a — )■ 0. 

The best fit point in a and tan /? is represented by 
the blue x and has a great overall fit of /v k, 0.5. 
Associated with the best fit, we find a small Mfj± that 
satisfies the current direct search bound fl9| . driven by 
the need to fit the large /i — !• 77 rate. In Fig. [21 we 
show the posterior mass distribution of the heavy states, 
which have the peak values different from the decoupling 
limit. Note that it is difficult to get such a large bb value 
of pvh{bb) = 1-3 indicated by the data. This is due to 
the already large branching fraction in the SM. Hence, 
no central value line is shown for that channel. However, 
this large value is dominated by the Tevatron Vh — > bb 
measurement. If the future rate measured at the LHC 
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FIG. 3: Posterior probability regions of la and 95% C.L. in 
the plane of 7^ and jd- Traditional 2HDM models (black cir- 
cles) the best fit of 2HDM-X (blue x) and 2HDM-D (middle 
black circle) are shown for comparison. 



is smaller as indicated by the CMS measurement, agree- 
ment within this channel for the 2HDM-X will improve. 
With a suppressed bb rate, another means to enhance 
the h "f-f rate is available as the total width decreases, 
thereby increasing the branching fraction to the other 
decay modes. However, we find this occurs within the 
2HDM-X infrequently. 

In Fig. [3l we show the posterior probability density in 
the plane of 'ju and 7^ . The most preferred values do not 
fit any of the traditional models. Of these models, the 
2IIDM-II has the most tension as it lies outside the 95% 
C.L. region. The lepton-specific (2HDM-L) case provides 
the best fit among the traditional models due to the in- 
dependence of the TT and qq couplings. This is verified 
by the reduced distribution by channel tabulated in 
Table |lTl We note that the values of 7„ = 7^ = 7r/4 ap- 
pear to give quite good fits for the democratic 2IIDM-D 
model. Among all the models we consider, the rr chan- 
nel provides the most tension. However, the 2HDM-L, 
2HDM-D, and the best-fit point of the 2HDM-X model 
provide excellent fits to the lepton data. 

As shown in Table HIl the best fit over all the data is 
given by the 2HDM-X, with the 2HDM-D not far behind. 
The pull from the data of the best fit point in the general 
model is indicated in Fig. S) We immediately see that the 
TT gives negligible rate while the 77 rate is well matched. 
It is worth mentioning that while the SM still gives an 
acceptable jv k, 1, the large tension within the rr and 
77 channels that persists even for the more traditional 
2HDM warrant consideration of the more general model 
as a potential explanation of the data. 
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FIG. 4: Channels measured at the LHC and Tevatron with 
best fit points from the 2HDM-X model denoted by x . 



TABLE II: Reduced x values for key channels in our anal- 
ysis. For each decay channel, we sum over the production 
mechanisms. 
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0.48 
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0.23 


0.59 
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0.59 
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0.20 
0.18 


0.23 
0.20 


0.83 
0.60 


1.31 
1.31 


0.53 
0.50 



Conclusions. The 2HDM-X framework provides a rich 
setting in which to analyze the LHC Higgs-like data in 
the context of simple extensions of the SM to allow for 
two electroweak Higgs doublets. As one might expect, 
the general 2HDM-X provides the best fit to the current 
data, with the democratic and lepton-specific models not 
far behind. A light charged Higgs is indicated at the best 
fit point in the 2HDM-X. Further experimental measure- 
ments of the neutral Higgs boson properties and searches 
for the charged Higgs boson will verify or falsify models 
within the 2HDM-X framework. 

We note here that during the completion of our paper, 
another paper appeared that also analyzes general two 
Higgs doublet models with vanishing tree-level FCNC's, 
with a different parametrization (20| . 

Acknowledgements. We thank Daniel Chung for 
useful discussions. VB, LE, and GS are supported by the 
U. S. Department of Energy under the contract DE-FG- 
02-95ER40896. YB is supported by start-up funds from 



the Univ. of Wisconsin, Madison. YB thanks SLAC for 
their warm hospitality. 



[1] 

[2] 
[3] 

[4] 

[5] 

[6] 

[7] 
[8] 

[9] 
[10] 



G. C. Bianco, P. M. Ferreira, L. Lavoura, M. N. Re- 
belo, M. Sher and J. P. Silva, Phys. Rept. 516, 1 (2012) 

LarXiv:1106.0034' [hep-ph]]^ 

N. Craig and S. Thomas. laFXiv:1207.4835l [hep-ph] . 

D. S. M. Alves, P. J. Fox and N. J. Weiner, 

.arXiv:1207.5499i [hep-ph]. 

S. Chang, S. K. Kang, J. -P. Lee, K. Y. Lee, S. C. Park 

and J. Song,'arXiv:1210.3439 [hep-ph]. 

R. D. Peccei and H. R. Quinn, Phys. Rev. Lett. 38, 1440 

(1977). 

V. Barger, H. E. Logan and G. Shauglmessy, Pliys. Rev. 
D 79, 115018 (2009) |arXiv:09 02.0170 [hep-ph]]. 
G. Bernardi and M. Herndon, arXiv:1210.0021 [hep-ex]. 
A. Rich and P. Tuzon, Phys. Rev. D 80, 091702 (2009) 
pxXiv:0908.1554 [hep-ph]]. 

J. L. Diaz-Cruz, A. Diaz-Furlong and J. H. Montes de 
Oca. larXiv: 1010 .0950 [hep-ph]. 

J. F. Gunion and H. E. Haber, Phys. Rev. D 67, 075019 



[11 

[12: 
[13 
[14 

[15: 
[16: 

[17 
[18 
[19 

[20: 



(2003) hep-ph/0207010 . 

S. Davidson and H. E. Haber, Phys. Rev. D 72, 
035004 (2005) [Erratum-ibid. D 72, 099902 (2005)] 

hep-ph/0504050 . 

I. Low, J. Lykken and G. Shaughnessy, larXiv: 1207. 10931 

[hep-ph]. 

G. Aad et al. [ATLAS C ollaboration], larXiv: 1207.03191 
[hep-ex] ; larXiv: 1207.72141 [hep-ex] . 

[CMS Collaboration], CMS PAG HIG-12-020; CMS PAG 

HIG-12-015; CMS PAG HIG-12-018. 

[Tevatron New Physics Higgs Working Group and CDF 

and DO Collaborations], arXiv: 1207.0449 [hep-ex]. 

M. Ciuchini, G. Degrassi, P. Gambino and G. F. Giudice, 

Nucl. Phys. B 527, 21 (1998) hep-ph/9710335: . 

J. Beringer et al. [Particle Data Group Collaboration], 

Phys. Rev. D 86, 010001 (2012). 

Y. Amhis et al. [Heavy Flavor Averaging Group Collab- 
oration], arXiv: 1207. 1158 [hep-ex]. 

[LEP Higgs Working Group for Higgs boson searches and 
ALEPH and DELPHI and L3 and OPAL Collaborations], 
hep-ex/0107031 

S. 



W. Altmannshofer 
larXiv: 1210.24651 [hep-ph] 



Gori and G. D. Kribs, 



